Физико-математические науки 








УДК 539.3 ООТ: 10.12737/1276 


Распределение контактных напряжений под круговой пластиной, лежащей 
на мягком слое” 


Б. И. Митрин, С. С. Волков 
(Донской государственный технический университет) 


Рассматривается построение эффективного аналитического решения осесимметричной контактной задачи о 
взаимодействии круглой пластины с двухслойным упругим полупространством. Для этого используется дву- 
сторонне асимптотический метод. Под действием нагрузки и реакции со стороны слоя пластина изгибается. 
Решение задачи получено для различных толщин слоя и значений параметра гибкости пластины. Рассмотре- 
ны случаи существенного отличия свойств между слоями полупространства. Для этих случаев построены 
аппроксимации трансформант ядра интегрального уравнения задачи высокой точности. Проведён контроль 
точности полученных результатов путём сравнения с известным решением задачи об изгибе круговой плиты 
на упругом слое, лежащем на недеформируемом основании. Исследовано влияние действия распределённой 
нагрузки на гибкие и жёсткие пластины, в зависимости от толщины слоя и его жёсткости по отношению к 
подстилающему полупространству. 
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Введение. Задача об изгибе пластины на упругом изотропном и однородном основании рассмат- 
ривалась в работах [1, 2]. Решение строилось путём представления контактных напряжений в ви- 
де степенного ряда, с последующим определением коэффициентов разложения из бесконечной 
алгебраической системы уравнений. 

Методом ортогональных многочленов такая задача решалась в работах [3, 4], а методом 
коллокации по чебышёвским узлам — в работах [5, 6]. При этом возникала необходимость 
построить решение некоторых бесконечных систем линейных алгебраических уравнений и 
ставилась проблема исследования сходимости полученного решения к точному. В работах [7, 8] 
для решения задачи применялись асимптотические методы типа «больших Л» и специальных 
ортогональных многочленов, что позволило получить основные характеристики решения в 
нескольких формах, каждая из которых эффективна в своей области изменения характерных 
параметров задачи. 

Отметим, что большинство известных решений эффективны только для жёстких пластин. 
И очень немногие, в частности, представленные в [7, 8], эффективны или для гибких, или для 
жёстких пластин. 

Интерес к решению задачи и её актуальность сохраняется и в настоящее время. Так, в ра- 
боте [9] решение строилось с использованием разложения напряжения в двойной ряд Фурье. 
Аналогичный подход использовался в работе [10]. Апагеа В. О. $Има с соавторами развил числен- 
ные методы решения задачи [11]. 

В настоящей работе развивается подход, основанный на двусторонне асимптотическом 
методе решения парных уравнений [12], позволяющий получить приближённое решение задачи в 
единой аналитической форме, применимой во всём диапазоне изменения геометрических и физи- 
ческих параметров задачи как для гибких, так и для жёстких пластин. 

Постановка задачи. Круглая пластина радиуса А и постоянной толщины р лежит на границе Г 
упругого полупространства © , состоящего из однородного мягкого слоя (покрытия) толщины Ни 


* Работа выполнена при поддержке Министерства образования и науки РФ (ГК № 11.519.11.3028, соглашение 
№ 14.132.21.1693), а также гранта РФФИ № 13-08-01435-а. 
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однородного полупространства (подложки). С полупространством связана цилиндрическая 
система координат /’, ф, 2'; ось 2’ нормальна плоскости Г и проходит через центр пластины. 


Под действием осесимметричной распределённой нагрузки р’(г’) и реакции со стороны слоя 


пластина изгибается (см. рис. 1). Функция ии” (^’) описывает прогибы плиты. 





Рис. 1. Схема для постановки задачи 


Введём обезразмеривающую замену переменных: 
л=Н/ Ву г=г/В; и (г) =и (ГВ, 
р” (г’)=р(г)рк?, 9* (г’)=9(г)9к?, 
где р — цилиндрическая жёсткость пластины, Па*м®. 
Коэффициенты Ламе двухслойного полупространства Л(2) и М(2) меняются по следую- 
щему закону: 


ме) [№ -1<2<0, 
Л, -ю<2<-1 (1) 
ме) ео, 

М, -®<2<-1. 


Модуль Юнга основания равен ЁЕ, =ВЕ, (-1). Безразмерный параметр В определяет скачок 
упругих свойств на границе «слой — подложка». Увеличение значения параметра В говорит об 


увеличении жёсткости подстилающего полупространства. 
Уравнение изгиба пластины имеет вид [14]: 


иг (г) = р” (г')-9° (г), 0<г’<В, 
где № — дифференциальный оператор изгиба пластины; р” (^’) — внешняя нагрузка, Па; 


49° (г’) — контактные напряжения под плитой, Па; и” (г’) — функция прогиба, м. 


Функция прогибов пластинки должна удовлетворять следующим граничным условиям: 
аи ху аи а 
а? га 


= 0, (Уи/) =0. (2) 
Условия (2) соответствуют свободному краю пластины, у — коэффициент Пуассона пластины. 








аг 


г=1 Г 
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При сделанных предположениях решение задачи сводится к решению следующей системы 
уравнений [13]: 


и (г) =р(г)-9(г), 0<г<1, (2) 
(в) (ел), (а) да = $и (г), 0<г<1, 

(3) 

[9() Л (аг)да=0, г>1, 

2 
где Ё, = [= 1%] ; [(а^) — трансформанта ядра. 
Параметр $ характеризует изгибную жёсткость пластины: 

5=ожр,, (4) 


К —1 
где © характеризует упругие свойства слоя на поверхности (© =2М, (Л, +М,) (Л, +2М,)_, Па). 
Парное интегральное уравнение (3) устанавливает связь между контактными напряжени- 
ями и прогибами пластины. Связь функции О (а) с контактными напряжениями 4(/) имеет вид: 


)- [9 ‚ (ар)рар, )- [0 (6) 2 (ог) ада. (5) 


Построение решения. Представим функцию прогибов в виде ряда по формам собственных 
колебаний круглой пластины со свободным краем, аналогично работе [14]: 


и(г)- Хи,» (^), „м р)рар, (6) 
где Ф, (г) = 2, ф„ (г) = А, [2 (К„г)- 2, (к) (К) 1 —. 


Постоянные А, и К,, для т =0,1,...,10 приведены в [14]. 


Учитывая линейность задачи, разложим функцию контактных напряжений 4(/г) в ряд сле- 
дующего вида: 


9(г)= Уи, (г), 0</<1. (7) 


Функции 9» (г) находятся из парного интегрального уравнения (3) с правой частью вида (6). 

Для построения решения парного интегрального уравнения (4) воспользуемся известным 
решением, полученным в работе [15], в случае, когда правая часть (4) представима в виде ряда 
по функциям Бесселя и учтём, что трансформанта ядра [ (а) может быть аппроксимирована сле- 
дующим выражением [12]: 

М 2 2 
Нея и) Поль" а„В, С. (8) 


Таким образом, решение парного интегрального уравнения (4) с правой частью (7) можно полу- 
чить в виде: 


о (Г) (0) п) + Усмнал" ) 


ди (г) =21"А„5 = (АК) (г, Жи) 7, (Ки) 1" (К) (РАК, Ч (г, К) + 


кУсич (ал Е ь, 


1=0 
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Постоянные С” определяются из системы линейных алгебраических уравнений: 
>. а(а,^",6,^") +1 (0) 6" =0, К=1,2,...,М 
УС" а(а,^", 6,^") +В (к А“) - Л (К) (К„)у (К, 9“) =0, 
р К=1,2,...,М; т=1,2,...., 
где а(а,Б) = [азва + 6 спа] (5? - а)" ‚ В(а,Ь) = 1; (ла)а(ва,Ь), у(а, Ь) = 1; (Ра)а(а,Ь). 
В свою очередь, контактные напряжения 4” (г) и "' | можно представить в виде ряда: 


9» (г) = т Ф, ( = [а р)рар, (9) 


р(г) = Хр» (г), = [РФ (р)рар. (10) 


Тогда, подставив разложения (7), (9), (10) в (3), имеем бесконечную систему линейных 
алгебраических уравнений для определения коэффициентов и’„, которую запишем в 
каноническом виде: 


ии, Е” ах! т =0,1,2,...; (11) 
750 
м 
= 21 *А,5 | (АК, Е - Л (к) (к) (РАК, ву + о, 


р = [Чоу (р)рар = А, [х (к, в») 2, (К) (к) (к, К], 


Ку = | (р,к,)Ф„ (р)рар = А, [Х (к, К) Л (К, (к„)х(ж, „|, 

ы [9 а,^")ф„ (р)рар = А, [Х (а,^7, к„)- 2 (к„) 1 (к, )х (ам, ж,)|. 
Здесь ре т = 0,1,2,..., 

(ав) 6 [совазть + 8 [ыт(а-6{а-Ь)"-чи(а+Ь)а+6) "| 


Х (а, Ь) = (а? + ь?) (азпасозр +р т спа). 
В частности, 


= 85 6 (о) АК» (п К» — Л, (к) (ке) (к) 5.157 | 
О] 
9 =21 "2,5 | [у (АК, Ку’ пк, — 2, (К, 1" (К, у (Ак, кли (К, ) + 


> с!а; эн(а,^- .| 


Используя метод редукции, сведём решение бесконечной системы алгебраических уравнений (11) 
к решению системы линейных алгебраических уравнений следующего вида: 


М 
ик +У и ЕТ =р:х =... М, (12) 
75 
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После определения коэффициентов и/„, (т = О. ле› М) для фиксированного значения М из (12), 


находим контактные напряжения 


М 
9(г)= У» и„9» (г), 0<г<1 (13) 
т=0 
и прогибы пластины 
М 
и (г) = \и„Ф„ (г), 0<г<1. (14) 
т=0 


Численные результаты. Рассмотрим случай мягкого упругого однородного слоя, лежащего на 
более жёстком упругом основании. Предполагаем, что параметр В принимает значения 2, 100, 
1000; коэффициент Пуассона слоя равен коэффициенту Пуассона подложки: \, =\, =0,3. 
Считаем, что слой жёстко сцеплен с основанием (аналогично можно рассмотреть случай, когда 
слой свободно лежит на подложке). 

На рис. 2 приведены графики трансформант интегрального уравнения (4) и их аппрокси- 
мации, задаваемые формулой (8), в случае отличия упругих свойств слоя и основания в 2, 100 и 
1000 раз. Причём рис. 2, а соответствует случаю отличия упругих свойств покрытия и основания в 
2 раза, 2, 6 — в 100 раз и 2, в— в 1000 раз. 


(и) 1 





0 2 4 цб $ 10 


Рис. 2. Трансформанты ядра интегрального уравнения и её аппроксимации 
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Для оценки погрешности аппроксимации трансформанты выражениями вида (8) использо- 
валась следующая формула: 

ы (и) (и) 
(и) 
Максимальная погрешность аппроксимации трансформанты ядра не превышает 3,4 % для 

графиков на рис. 3. 

В случае, когда подложка в 2 раза более жёсткая, чем покрытие, погрешности аппро- 
ксимации менее 3,5 % можно добиться при Л = 1 (соответственно 1,4 % при № = 8), а в случае, 
когда основание в 1000 раз более жёсткое, чем покрытие, для снижения погрешности до 3,1 % 
понадобилось ЛМ = 20. 

В работе [16] показано, что по мере роста значений параметра В трансформанты ядер 


А, (и) = 100%. (15) 








приближаются к предельному графику, соответствующему случаю недеформируемого основания. 
Ниже покажем, что значения контактных напряжений для В=1000 и недеформируемого 
основания [17] близки друг к другу. 

Рассмотрим случай, когда на пластинку действует равномерно распределённая нагрузка. 
Упругий однородный слой мягче основания в 1000 раз. 


Ч Ч 

0.90 

р) 
т 0.85 

0.90 
0.80 

0.88 
0.86 0.75 
0.84 0.70 





0 0.2 0.4 0.6 0.8 Г 0 0.2 0.4 0.6 0.8 г 


Рис. 3. Контактные напряжения для однородного слоя, лежащего на упругом (В = 1000 , сплошная линия) 
или недеформируемом (прерывистая линия) основании 


На рис. 3 сплошной линией изображены графики контактных напряжений, построенных 
для аппроксимации трансформанты ядра при В =1000, график которой построен на рис. 2, в. 


Прерывистой линией обозначены результаты из работы [17], полученные для средней зоны 
значений параметра Л, как для гибких, так и для жёстких пластин (параметр $ = 0,1 и $= 4 
соответственно). Сравнивая данные результаты с полученными в этой работе, можно обнаружить, 
что их значения близки друг к другу, разница между ними не превышает 6,5 %. Заметим, что 
графики на рис. 3 построены двухсторонне-асимптотическим методом, для которого доказана 
сходимость метода как для больших, так и для малых значений характерного геометрического 
параметра задачи Л [18]. Анализ графиков на рис. 3, показывает, что метод охватывает и зону 
средних Л. Таким образом, можно говорить, что при построении достаточно точной 
аппроксимации трансформанты ядра интегрального уравнения (4), используя двусторонне- 
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асимптотический метод, можно получить решение задачи во всём диапазоне значений 
характерного геометрического параметра Л, как для гибких, так и для жёстких пластин. 
Рассмотрим случай В=100. Трансформанта ядра и её аппроксимация для закона, 


описанного выше, построена на рис. 2, 6. 

На рис. 4 изображены контактные напряжения для гибкой пластины ($ = 0,1). Как видно 
из графиков, для больших и средних значений параметра Л (диапазон изменения Л от 1 до 4), 
распределение контактных напряжений близко к значениям, соответствующим однородному по- 
лупространству. Для малых Л (Л =0,1; 0,25 ) графики контактных напряжений сильно отличаются 
по поведению и по значениям от соответствующих однородному полупространству. Таким обра- 
зом, можно сказать, что относительная толщина мягкого упругого слоя существенно влияет на 
распределения контактных напряжений под пластиной. 





9(г) 9(г) 
1.02 
0.97 
1.00 
ее 0.98 
| 0.96 
0.95 0.94 
0.92 
0.94 0.90 
0 02 04 06 08 0.95" 0 02 04 06 080.95 
а(г) ^А=0.1 (г) =0.25 
1. 1.5 
ты 1.4 
| 1.3 
1.2 1.2 
1.1 1.1 
т 1.0 
0.9 
0.9 0.8 
0.8 0.7 
0.6 
0.7 ь ь 
о 02 04 06 08 0.95 0 02 04 06 058095 


А=1 ^л=4 
Рис. 4. Контактные напряжения. Закон 1, $ = 0,1, № = 10, М= 10 
На рис. 5 изображены контактные напряжения для жёсткой пластины ($5 = 3). Анализ гра- 


фиков на рис. 5 аналогичен рис. 4. Отметим интересный качественный результат, заключающий- 
ся в том, что при малых Л (Л=0,1;0,25) максимум контактных напряжений в диапазоне 


ГЕ [0..0,95 ] достигается в центре пластины. 
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9 —_—_—_ а(г) 
0.98- 
| 0.98 
0.97- 
| 0.96 
0.96- 
1 0.94 
0.95- 
1 0.92 
0.94- № 
] 0.90 
0.93- \ / 
| т т т т 0.88 т т > г 
0 0.2 0.4 0.6 05 0.95" 0 0.2 0.4 0.6 0.8 0.95 
^=0.1 ^=0.25 
9(г) 
1.4- 
ЕЯ 
1.21 
1.11 
1.0- 
0.9- и 
0.8- р 
0.7- —“ 
0 02 04 06 080.95" 0 0.2 04 0.6 0.8 0.95" 
^=1 ^А=4 


Рис. 5. Контактные напряжения. Закон 1, $ = 3, М = 10, М= 10 


Заключение. С помощью двусторонне-асимптотического метода получено эффективное 
аналитическое решение осесимметричной контактной задачи о взаимодействии круглой пластины 
с двухслойным упругим полупространством. Найдены в аналитическом виде контактные 
напряжения и функция прогибов пластины. Построенное решение является двусторонне- 
асимптотически точным, то есть оно эффективно как в зоне больших, так и малых значений 
характерного геометрического параметра задачи (отношение толщины слоя к радиусу пластины) 
для гибких и жёстких пластин. Численные результаты, приведённые в работе, показывают, что 
полученные аппроксимации высокой точности для трансформанты ядра интегрального уравнения 
задачи позволяют получить решение, которое эффективно и в зоне средних значений 
характерных геометрических параметров задачи. 
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(Роп За Тесйптка! Упмег®Ку) 
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